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Abstract 

We consider the 3-body problem in relativistic lineal (ie 1+1 dimensional) gravity 
and obtain an exact expression for its Hamiltonian and equations of motion. While 
general-relativistic effects yield more tightly-bound orbits of higher frequency compared 
to their non-relativistic counterparts, as energy increases we hnd in the equal-mass case 
no evidence for either global chaos or a breakdown from regular to chaotic motion, 
despite the high degree of non-linearity in the system. We find numerical evidence 
for mild chaos and a countably infinite class of non-chaotic orbits, yielding a fractal 
structure in the outer regions of the Poincare plot. 
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Self-gravitating iV-body systems have long been of key importance in studying both 
stellar dynamics (small N > 2) and galactic evolution (large N). Solutions to their problem 
of motion, complicated for N > 2 even in nonrelativistic mechanics due to singularities and 
evaporation, are essentially intractable for iV > 1 in general relativity because of gravita- 
tional radiation. These difficulties are avoided in one-dimensional (or lineal) models of such 
systems, which describe a system of N parallel mass sheets. Studied by astronomers and 
physicists for several decades [1], they merit special attention [2, 4], since they not only 
admit a significantly greater level of computational and analytic analysis but can also be 
mapped to systems subject to experimental test [3]. 

We report here on the results of an investigation of the 3-body problem for a relativistic 
self-gravitating lineal system. Three is the smallest number of particles for which chaotic 
motions are possible and so such systems are of particular interest. Non-relativistically this 
system is equivalent to that of a single particle moving in a hexagonal funnel [5] or to that 
of a uniformly accelerated particle undergoing elastic collisions with a wedge [6]. While 
its dynamics are well-documented, ours is the first non-perturbative relativistic treatment 
of this problem. From the canonical formulation of lineal gravity minimally coupled to 
3 particles, we obtain an exact expression for the Hamiltonian (valid to all orders in the 
gravitational coupling k = 87rG/c 4 ), transcendentally expressed as a function of the four 
independent degrees of freedom. From this we are able to solve for the motion of the system 
under arbitrary initial conditions and study quasi-periodic motion and chaos in a general- 
relativistic self-gravitating system. 

For our lineal self-gravitating system we minimally couple iV point masses to gravity in 
2 spacetime dimensions 

(1) 

where g^ u and g are the metric and its determinant, R is the Ricci scalar, r a is the proper 
time of a-th particle, and we incorporate a scalar (dilaton) field \I/ since the Einstein action 
is a topological invariant in 2 spacetime dimensions. The action (1) describes a generally 
covariant self-gravitating system (without collisional terms, so that the bodies pass through 
each other) that is a generalization of Jackiw-Teitelboim lineal gravity [7], in which the 
scalar curvature is equated to a (cosmological) constant A; its equations of motion are 

iv„*V„* - ft ,„ f jV A Wy» - V 2 *) - V„V„* = kT? + ^J„„ (3) 
where the stress-energy due to the point masses is 
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and is conserved. Eq. (2) is a closed system of A" + 1 equations for which one can solve for 
the single metric degree of freedom and the A" degrees of freedom of the point masses. The 
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evolution of the dilaton field is governed by the evolution of the point-masses via (3). The 
left-hand side of (3) is divergenceless (consistent with the conservation of T^ u ), yielding only 
one independent equation to determine the single degree of freedom of the dilaton. 

For general N the canonical formulation of (1) has been studied previously [8] and a 
number of exact solutions to the 2-body problem have been obtained [9]. Writing the 

metric as ds 2 = — (N Q dt) 2 + 7 (dx + ^dt^j (with f dx^j describing proper spatial distance 

at fixed t), the action can be canonically reduced to I = f d 2 x |X)^Li Pai a S(x — z a ) — 

where the reduced Hamiltonian is H = J dxH = —\ j dx^fr" , with the overdot and prime 
denoting and ^ respectively. The field \I> = ^(x,z a ,p a ) is understood to be determined 
from the constraint equations, which become 



a 
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(5) 
(6) 



where n = x' and p a are the respective momenta conjugate to 7 and z a . We have chosen the 
momentum conjugate to \I/ to vanish and 7 = 1, fixing the frame of the physical space-time 
coordinates in a manner analogous to the (3 + 1) -dimensional case [10]. 

Remarkably eqs. (5,6) can be exactly solved when N = 3. The solution for ^ is a 
sum of exponentials, and for x a linear function of x in between each particle. Requiring 
finiteness of the Hamiltonian at x — ±00 and a consistent match of the solutions across the 
delta-functions at the particle boundaries yields 
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M t = + U = H - Mi- '(Y.,!^,). L* = (1 - 
Ilj<fc^i s ij s ik)Mi + Li, Zij = (zi — Zj) and Sij —sga(zij). For simplicity we have set A = 0. The 
discrete constant of integration e = ±1 flips sign under time- reversal; it provides a measure 
of the flow of time of the gravitational field relative to the particle momenta. 

Eq. (7) is an exact result and implicitly determines the Hamiltonian H as a function 
of the two independent coordinate degrees of freedom (p, A) and their conjugate momenta, 
which can be written as 



z l — z 2 
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- 2z 3 = V6\ 
2p 3 = V6p x 



(8) 
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choosing the center of momentum to vanish. The usual relations z a = and p a = — f^- 
yield the equations of motion, and it is straightforward to show that the Hamiltonian is 
time-independent when these equations hold. 
A post-Newtonian expansion [8] of (7) yields 
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in the equal mass case. The first term is the total rest mass and the next two terms are 
equivalent to the hexagonal- well Hamiltonian of a single particle [5], which we refer to as 
the hex-particle. We see that to leading order in c -2 the hexagon potential is modified to 
have parabolic edges, as well as a momentum-dependent steepening of its walls. For unequal 
masses the hexagonal symmetry becomes distorted, with two opposite corners moving in- 
ward, changing the relative arc-lengths of the sides. The full potential (obtained from (7) by 
setting p p = p\ = 0) has the shape of a hexagonal carafe with convex edges. 

The system's interesting dynamics are due to the non-smoothness of the potential along 
the lines p = 0, p — \f?>\ = 0, and p + \/3\ = 0, respectively corresponding to the crossings 
of particles 1 and 2, 2 and 3, or 1 and 3. We distinguish two distinct types of motion [6]: A- 
motion, corresponding to the same pair of particles crossing twice in a row (the hex-particle 
crossing a single sextant boundary twice in succession), and 5-motion, in which one particle 
crosses each of its compatriots in succession (the hex-particle crossing two successive sextant 
boundaries). We can characterize a given motion by a sequence of letters A and B (called a 
symbol sequence), with a finite exponent n denoting n-repeats and an overbar denoting an 
infinite repeated sequence. 

We numerically solve the equations of motion that follow from (7) in both the exact 
relativistic (R) system and its non-relativistic (N) limit (ie c — > oo). We impose absolute 
and relative error tolerances of 1CT 8 for the numerical ODE solvers. This yields numerically 
stable solutions for the values of H < Qmc 2 that we study; beyond this the ODE solving 
algorithm is unstable. We test stability by checking that the energy remains a constant of 
the motion to within a value no larger than 10~ 6 throughout. 

For both N and R systems we find three principal classes of motion. One class, with 
sequence B, describes annulus-shaped orbits encircling the origin in the p — A plane. In the 
second (pretzel) class the hex-particle essentially oscillates back and forth about one of the 
three bisectors in a combination of A and B motions. Chaotic orbits constitute the third 
class, where the hex-particle wanders between A- and 5-motions in an apparently irregular 
fashion. Such orbits eventually wander into all allowed areas of the p — A plane. 

Although the orbits of the N and R systems realize the same symbol sequences, important 
qualitative differences exist between them. For a meaningful comparison between each sys- 
tem we set the total energy H = Et and the initial values of (p, X,p p ) to be the same, setting 
the initial value of p\ to satisfy (7) in each case. As the parameter r\ = H/ (3mc 2 ) — 1 grows, 
we find that R trajectories have higher frequencies and extend over a smaller region of the 
(p, A) plane than their N counterparts. Choosing initial conditions so that A-motion takes 
place at A = we find that R trajectory patterns narrow with increasing r\ in the small-A 
region, whereas those for the N system do not. Figure (1) illustrates a typical example with 
7] = 0.5 for the pretzel case, in which the trajectories for the 3 bodies and their corresponding 
hex-particle are shown. For each, two of the bodies form a low-amplitude/high-frequency 
bound state that in turn oscillates at higher amplitude and lower frequency with the third. 
The 3-body N oscillations are parabolic in shape whereas there is shoulder-like distortion in 
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Figure 1: Regular pretzel orbits (run for 120 time steps) of the R and N systems with the 
corresponding 3-particle time evolution (truncated at 50 time steps). The collision sequences 
differ for these initial conditions and correspond to AB 3 (R) and A 2 B 3 (N). 

the R system due to the momentum-dependence of the potential (a feature seen previously 
in the 2-body case [9]) illustrating that under appropriate initial conditions two bodies can 
tightly and stably bind together in both the N and R systems. The hex-trajectories differ 
substantively, with the R pattern having 50% higher frequency, and extending over a slightly 
smaller range than its N counterpart. The R- "boomerang" also has an indentation on the 
left at A = that is absent for all corresponding N system boomerangs (not shown). 

We probe the global structure of the phase-space using Poincare plots, which we construct 
by recording the radial p r and angular momentum j each time the hex-particle crosses a 
hexagonal bisector, plotting their values as a series of dots in the {p r ,j 2 ) plane. Periodic 
trajectories appear as a finite series of dots, quasi-periodic trajectories as closed loops and 
chaotic ones as densely filled regions. Our numerical results (see fig. 2 ) for the N system 
agree with previous studies [5, 6, 11], which established that a region of phase space exists 
in which the hex-particle circulates regularly around the origin, bounded by a thin region 
of chaos outside of which there is a regular structure of periodic and quasi-periodic orbits. 
Remarkably, we do not observe a breakdown from regular to chaotic motion as f] increases, 
despite the high degree of non-linearity in the R system. The lower regions of the Poincare 
map clearly display the same pattern of series of circles as occurs in the non-relativistic 
case. However the Poincare plot for the R system is no longer symmetric with respect to 
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Pohcare Plot of (he N System 




Figure 2: The Poincare plot of the N system. The squares denote the parts of the plot 
magnified in the insets. 

p r = 0, but instead is 'squashed' towards the right-hand side (compare fig. 2 to fig. 3). 
This deformation is reminiscent of the situation for two particles, in which the gravitational 
coupling to the kinetic-energy of the particles causes a distortion of the trajectory from an 
otherwise symmetric pattern [9]. Since (7) is invariant only under the discrete symmetry 

e) — > (— Pi, — e) (and not under p^ — > —pi), we obtain for the choice e = ±1 a distortion 
towards the lower right/left-hand side of the plot relative to its N counterpart. 

No sizeable connected areas of chaos are present in either the R or N systems, indicating 
that most trajectories are effectively restricted to move on two-dimensional surfaces in phase 
space. We find for all r\ within the range investigated that chaotic orbits exist separating the 
annulus and pretzel regions. We therefore conjecture that the R system is not integrable, 
since its N counterpart is known to be non-integrable [6] . Nonetheless, clearly some underly- 
ing feature enforces considerable structure on the phase space, preventing KAM-breakdown 
to global chaos typical of such Hamiltonian systems [12]. 

In general pretzel-type orbits display a remarkable richness of dynamics. As the ini- 
tial angular momentum of the trajectory in question decreases, the number of successive 
A collisions increases before the hex-particle sweeps around the origin in a -B 3 -sequence, 
corresponding to a 180-degree swing of the hex-particle around the origin . For example 
AB 3 (the simplest sequence after B) corresponds to a boomerang-type orbit and appears 
as two circles on the Poincare section. The next simplest sequence is A 2 B 3 , which corre- 
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Figure 3: The Poincare plot of the R system at rj = 0.75 and e = +1. The upper right inset 
provides a close-up of the chaotic region at the top of the diagram; it is now considerably 
narrower than for lower values of rj. The lower-right inset is a close-up of the structure in a 
pretzel region in the lower right of the diagram. 
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sponds to a bow-tie like orbit, and generates three slightly smaller circles on the Poincare 
section. For each of these patterns a family of orbits exists corresponding to different widths 
of the 'bands' of phase space that the trajectory covers, and correspondingly different radii 
of circles in the Poincare section. Between these regions, the orbits' sequences are mixtures 
of AB 3 and A 2 B 3 . This reasoning can be extended to more general A n .B 3m -motions. We 
conjecture that the only allowed non-chaotic orbits - relativistic and non-relativistic - are of 
the form Hnrnj {A n B 3m ) 3 with n,m finite, corresponding to increasingly complex weaving 
patterns. 

If the set of integers lj is finite, then the sequence is regular, leaving bands of phase 
space untravelled, and appearing as a series of closed crescents or ellipsoids on the Poincare 
section. If, however, the sequence of integers lj never repeats itself, then the trajectory will 
fill the available phase space densely, appearing as a wavy line on the surface of section. We 
conjecture that there is a 1 — 1 correspondence between rational numbers and periodic orbits 
in this region of phase space, both non-relativistically and relativistically. This would give 
the lower section of the Poincare plot a fractal structure as the patterns of circles, ellipses 
and lines is repeated on arbitrarily small scales as the hex-particle's angular momentum 
approaches zero. 

The 3-body relativistic system provides a new theoretical laboratory for studying the 
interplay between relativity and chaos. While the dynamical structure of the relativisitc 
system bears some resemblance to well-known features of its non-relativistic counterpart 
(such as the self-similar structure of the pretzel class and the lack of energy dependence of 
the orbit topology), significant differences arise in the shape and frequency of the orbits. A 
full study of the large-r/ regime should provide us with interesting new information to this 
end. 

This work was supported by the Natural Sciences and Engineering Research Council of 
Canada. 
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